Prepared for submission to JCAP 



Transporting non-Gaussianity from sub 
to super-horizon scales 

David J. Mulryne 

Astronomy Unit, School of Physics and Astronomy, Queen Mary University of London, 
London, El 4NS, UK 

E-maih d.muhyne@qmul.ac.uk 

Abstract. We extend the 'moment transport method' for calculating the statistics of infla- 
tionary perturbations to the quantum phase of evolution on sub- horizon scales. The quantum 
transport equations form a set of coupled ordinary differential equations for the evolution of 
quantum correlation functions during inflation, which are valid on sub- and super-horizon 
scales, and reduce to the known classical transport equations after horizon crossing. The 
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1 Introduction 



The non-Gaussian statistics of the curvature perturbation, produced during inflation, have 
become a key test of inflationary models. Methods to calculate these statistics are now well 
established. Typically the calculation is split into two parts, one for the quantum 'sub- 
horizon' regime, and one for the classical 'super- horizon' regime. It is common to perform 
the former part of the calculation using standard techniques of quantum field theory. First 
the interaction picture is employed, then initial conditions are provided by assuming linear 
Fourier modes whose wavelength is far smaller than the cosmological horizon are in the 
Bunch-Davis vacuum state, and finally same-time correlation functions are calculated in 
the interacting vacuum. This calculation is often described as the In-In formalism applied 
to inflationary cosmology (see for example Refs. [1-5] for calculations of the inflationary 
bispectrum). It is usual to evaluate these quantum correlation functions shortly after the 
Fourier modes of interest cross the cosmological horizon, where for single field models (, and 
hence the correlation functions of C, become constant [6, 7]. 

In multiple-field models, a complication is that C can continue to evolve. In principle 
there is no barrier to evaluating the correlation functions at later times using the In-In 
formalism, thus accounting for this evolution. But in practice this becomes extremely difficult 
analytically (see Ref. [8] for recent progress in that direction), and moreover the integrals 
which one is required to evaluate are not particularly well suited to numerical integration. The 
usual procedure, therefore, is to truncate the full quantum calculation at horizon crossing, to 
interpret the correlation functions as classical statistical quantities, and to employ classical 
methods to follow the statistics from that point on [2, 9]. This is valid since the dynamics is 
expected to be very close to classical in this regime. 

Despite the fairly mature state of existing techniques, it is still useful to explore new 
methods to calculate the statistics of inflationary perturbations for several reasons. First, 
with the imminent expected influx of data, computationally efficient algorithms are urgently 
needed to precisely compare the predictions of inflationary models with observations. Sec- 
ondly, in this connection, it would be useful to have a computationally efficient method for 
which a splitting into sub- and super-horizon regimes was not required. This will be invalu- 
able for many models, such as those in which the inflationary trajectory bends during horizon 
crossing, where such a splitting will not yield accurate results. Finally, there is considerable 
interest in understanding the exact relation between sub-horizon and super-horizon methods. 
For example, in order to calculate the error incurred by using purely the latter as is com- 
monly done in studies of multiple field inflation. One approach to understanding this link 
is to consider how classical equations emerge from the quantum system via the dynamical 
renormalisation group flow [8]. In this paper we pursue a complementary route. 

The most widely applied super-horizon approach is the 6N formalism which utilises 
the separate universe approximation [9-13]. Recently another technique called the 'moment 
transport method' has also been extensively explored [5, 14-16] (also see [17-19] for related 
earlier work). One can show that objects which provide formal solutions to the transport 
system (see §3) are equivalent to the coefficients of a 5N style Taylor expansion [20, 21], 
and hence the methods are formally equivalent. The transport approach, however, has the 
computational advantage of providing a simple set of coupled evolution equations for the 
classical two-point, three-point, and higher, correlation functions, or for the Taylor expansion 
coefficients [21]. Such equations provide important analytic insights into the evolution of the 
statistics [20, 22, 23], and are amenable to a straightforward numerical implementation (see. 



- 2 - 



for example, Refs. [24, 25] for studies utilising numerical implementations of the transport 
equations) . 

In this paper, our aim is to show how the transport equations of motion can be extended 
to sub-horizon scales to provide a set of coupled ordinary differential equations (ODEs) for 
quantum correlation functions^, which are valid in both sub- and super-horizon regimes, 
and tend to the classical equations after horizon crossing. This offers an alternative route to 
calculating the statistics of C- Moreover, since differential equation formulations of a problem 
are often easier to implement numerically than integral formulations, the quantum transport 
equations will likely also be of benefit for numerical implementations of the full quantum and 
classical system. Though we defer such a implementation to future work [27]. 

The layout of the present work is as follows. In §2 we review the transport approach 
to classical statistics, and discuss how the transport equations look for equations of motion 
in which spatial gradient terms have not been discarded, as they typically are on super- 
horizon scales where transport methods have previously been applied. In §3 we review some 
analytical properties of the moment transport equations and related transport equations. 
Then in §4 we show how the formalism can be further extended to regimes in which the 
field perturbations behave as quantum objects, as they do on sub-horizon scales. We also 
show how formal integral solutions to these quantum equations are connected to the integral 
expressions of the In-In formalism. In §5, we explore these results using an explicit example 
of canonical single field inflation. Finally, in §7, we show that the transport expressions for 
the two-point function in multiple field settings recover known methods originally described 
by Salopeck et al. [28]. We conclude in §8 and discuss future directions. 

2 Transporting classical statistics 

The simplest way to derive the transport equations is to begin with the evolution equations for 
cosmological perturbations themselves^. In Fourier space, we can utilise a form of DeWitt 
notation to write these perturbation equations in a highly compact form [20, 21]. In this 
notation, x^' represents the Fourier modes of all fields and field velocity perturbations defined 
on flat hypersurfaces and evaluated at some time . A primed index, for example a', on 
any object indicates a field/field velocity label a and a momentum label k(x, such that for 
example 

Xa' = a;«(ka) . (2.1) 

a therefore runs from 1 to 2n where n is the total number of fields, and for odd a, x^ii = 
5(j)[a+i)' /2-, while for even a, Xa' = (5</'a'/2- Repeated primed indices imply that the fields 
and field velocities are summed over and an integration is performed over the suppressed 
Fourier space dependence, with measure d3fca/(27r)3. Repeated non-primed indices indicate 
just summation over the field label, leaving the Fourier mode unchanged as the sum runs 
through each value of a. Below we will see an explicit example of this. 

^Similar equations are used in other settings (see [26] and references therein), and provided inspiration for 
looking at this approach in the inflationary setting. 

■^One can also begin with the partial differential equation for the joint probability distribution of inflationary 
perturbations and decompose this equation into ordinary differential equations for each moment or cumulant 
[14, 15]. 

^Or any other complete set of perturbations which completely parametrise the perturbed spacetime, for 
example the curvature and isocurvature modes, can be used. In this paper we will always consider fields and 
field velocity perturbations. 
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In this notation, perturbations obey the equations of motion 
dx / 1 

— ^ = Ua'P'Xpi + -Ua'P'-^i {xp'Xy - (x/3'Xy)) , (2.2) 

where the u coefficients can be derived from cosmological perturbation theory (for a recent re- 
view see Ref. [29]), or on super-horizon scales from alternative methods utilising the separate 
universe approach [14, 15, 17, 18, 21]. The form ofua'/3' is 

Ua'i3' = (2vr)^'Ua/3(A;„)(5(k« - kp) , (2.3) 

where we sometimes describe Ua/sika) (distinguished by the unprimed indices) as the kernel 
of Ua'/3', and, for example, 

Ua'p'Xp' = j d^/cpnQ,/3(A:a)(^(kcc - k|3)x/3(kp) , 

= Ual3{ko,)xp{\<ia) , (2.4) 

where on the right hand side of the final expression the only sum is over the repeated index 
/3, and kg^ is simply one particular Fourier mode, which we could have equally well labelled 
A;, /ci, A;2, k' or k" etc. as is common in the literature. Likewise, the form of Ua'p>'y' is 

Wa'/3'7' = (2vr)^Ua/3T,(A;a, A;p, ky)5{\<ia - k|3 - k^) , (2.5) 

and for example, 

Ua'p'^'Xpix^f = -j-^—^ I I d^%d^A;yMQ/3^(A;<x, A;(3,A;Y)5(ka - kp - kY)x/3(k|3)j;^(kY) , 



(2vr)= 
1 

(2^ 



d^fc|3tia;3^(A;„, A;(3, A;a - A;p)2;/3(k|3)x^(ka - kp) . (2.6) 



The form of u coefficients given in Eqs. (5.1) and (5.2) are such that after substituting into 
Eq. (2.2), the various terms take the expected form, which follows from Fourier transforming 
the original real space equation of motion. In particular the second term on the RHS of 
Eq. (2.2) involves a convolution over all scales. In general the kernels of the u coefficients are 
functions of background quantities and can also carry Fourier mode dependence as indicated 
above. For super-horizon regimes they become independent of any scale dependence. In §5 
we will see explicit examples of the u coefficients which arise from cosmological perturbation 
theory. 

The notation involving primed indices we have introduced in this section may seem un- 
necessary. However, we can already see from Eq. (2.2) that it allows complicated expressions, 
involving for example convolutions, to be written in a compact form utilising an extended 
summation convention. This will become increasingly advantageous for manipulations as we 
proceed. 

2.1 Evolution of statistics 

In order to turn the evolution equations for field fluctuations into evolution equations for 
their statistics, we employ the simple principle that expectation values obey the equation 



dt \ dt 
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where A is any function of perturbations which has no exphcit time dependence. As discussed 
in more detail in Refs. [20, 21], this foUows by considering the fact that an evolving classical 
distribution of perturbations must conserve probability. 

By defining the two and three-point correlations functions as 

Sq,'/3' = {Xa'Xi3') , (2.8) 

and 

aa'13'y = (xa'X^'Xy) , (2-9) 

respectively, identifying (A) with Sq'^' and then with Ua'/s'y, and finally utilising Eqs. (2.7) 
and (2.2) and the chain rule, one finds a set of coupled equations 

dSo,/^' 

— = Uq^i^iYj^i pi + Upi-fiYi^i a' + . . . , 

^"^^^ ^ = Ua'X'axipi^i + Uq'a'/x'Sa'/S'S/x'^' + (a' /3' 7') + . . . , (2.10) 

where the dots indicate we have truncated the equations to leading order, and the arrow that 
there are two additional terms formed by cyclic permutations. Here we have considered only 
the two and three-point correlation functions, and for simplicity we will restrict ourselves to 
these throughout this paper. Similar expressions exist for the four-point function [21] and 
can be calculated for the higher n-point functions in the same manner. These equations allow 
known classical statistics at some point in time to be propagated forward to a later point 
at which we wish to evaluate them (for example between horizon crossing and the end of 
inflation) . 

2.2 Equations for correlations are simpler than equations for perturbations 

The above equations are complicated because they involve integrations over k space, which 
follow from the convolutions which naturally appear in Fourier space perturbation equations 
beyond linear order. Correlation functions in Fourier space, however, are always accompanied 
with delta functions, and at any given time, t, we have 

Ya'P' = {2^fT.^p{kc,)5{\^o, + kp) , (2.11) 

and 

aa'/3'7' = {2Tifaap^{ka, kfi,,ky)5{\<io^ + kp + k-y) , (2.12) 

where we refer to '^aiii^cx) and cta/B^/iko^, k^, ky^ as the kernel of the two- and three-point 
functions respectively"^. We find that 

^ = Ma^(/Ccc)^7/3(fea) + Uf}j{kc()TijQ,(k(x) + . . . , 



dt 

daQ,/3^(fca, fcp,fcy) 
dt 



= Uax{koc)axi3j{ka, k^,ky) + UaXfiika, k^,ky)^xp{k^)T.^^{ky) 

+ (a^/3^7) + ... , (2.13) 



where once again we give only the leading order terms in each equation. In the final expression 
it should be understood that moving the position of free indices when performing the cyclic 

*Note that for convenience, the notation is slightly different to Ref. [21]. 



- 5 - 



permutations, a — /3 — )• 7, also moves the positions of the associated Fourier modes fca — )■ 
— >• ky- An exphcit example of a system of this type is given in §5. These equations 
exhibit the feature that one does not need to perform a convolution over all scales in order 
to calculate the evolution of aaP'yikcc, kp,, ky), as one does in Eq. (2.2) for the evolution of 
XaO^a)- In this sense Eqs. (2.13) are simpler than Eq. (2.2). 



3 Other properties of the transport system 

The transport system, given by Eqs. (2.10), has many attractive properties. In particular, 
one can introduce objects, closely related to the correlations themselves, which allow formal 
analytic solutions of the transport hierarchy, using path-ordered exponentials (see Refs. [20, 
21]). Here we simply quote the results, and point out how they appear in a scale dependent 
setting. 

The solutions of Eqs. (2.10) from some initial time to to some later time t are given by 

= ra/j/r^/j/Sj/j/ , (3.1) 

+ (a'^/3'^7')] , (3.2) 

where we have adopted further compact notation according to which a Greek subscript indi- 
cates evaluation at the later time t as before, and a Roman index indicates evaluation at the 
initial time to, and where we have presented only the leading order terms in these solutions. 
Mixed index objects have dependence on both times. Here we have introduced a new function 
Tcj/j/ which is the propagator of the system, r^/j/ satisfies its own transport equation 

dTa'j' 

= Ua'p'i- p'i' , (,3.3j 
and we have also introduced Ta'i'j' which obeys the transport equation 

d^a'i'j' 

^ = Ua'l3''^l3'i'j' + Wa'/3'7'r/3'i'r7'j' • (3.4) 

The initial conditions are TJ^/*, = 5a'i' and T™/*,^., = 0, which follow from their definition in 
Eqs. (3.1) and (3.2). These objects can be solved in terms of a path ordered exponential to 
find 

Ta'i' = "Pexp dt' Ua'f{t')^ Syi' , 

^a'i'j' = / dt'Ta'fi'Ufi/c,'^'{t')Ti;,i'T^'ji , (3.5) 
J to 

where the indices marked with a tilde, indicate they are associated with the intermediate 
time in the integral t' . 

In terms of these solutions, Eq. (3.2) becomes 

"a'/3'7' = ^a'i'^iS'j'^y'k'Oii'j'k' + / dt' [Ta'fi'Ufiip'5-iI!,fS/u'T,Y^i + (ck' —)•/?'—)• 7')] , (3.6) 

Jto 

where we have introduced mixed indexed two-point correlations, Sa'j' = Ta'k''^k'i', which by 
construction satisfy the same equation of motion as T^'i'- 
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Once again the expressions above can be simplified, since consistency implies that 



r^,i, = {27^yTa^{K)6{K - ki) , (3.7) 

and 

Ta'i'j' = {2TrfTaijikcc, k, kj)5(ko, - ki - kj) . (3.8) 

As in the previous sections, therefore, all integrations over Fourier modes can be explicitly 
performed in the equations of motion (3.3) and (3.4), and in the integral solutions. We find, 
for example, 

^ai{ka) = 7^ exp dt' nQ^(fca,i')^ '^7* > 

^aijika, ki, kj) = / dtTap,{koi)ufiDaika,ki,kj)Ti^i(ki)Taj (fcj), (3.9) 
Jto 

and 

"a/37 = 'raiikaWl3j{kp)T^kiky)aijkika,kp,ky) 

olJU fiver (koc, k^,ky)T.i3i^{k^)T.j^{ky) 

Jto 

+ (a^/3^7)]. (3.10) 
One can verify that direct differentiation of Eq. (3.10) gives the evolution equation (2.13). 

4 Quantum transportation 

We now turn to quantum transport equations. Our aim is derive analogous equations and 
solutions to those above, but valid even when the quantum properties of the fluctuations are 
important. This will mean that if values of inflationary correlation functions are known at 
some initial time, even if that time is long before horizon crossing, they can be propagated 
forward to some later time, such as end of inflation. 

To do this we first promote dcpaO^a) to an operator dcpaO^a)- This will be accompanied 
by a canonical momentum = dC/d6<j)a, where C is the Lagrangian which leads to Eq. (2.2). 
5(j)a and pa form a conjugate pair obeying the usual commutation relations [Scpa'^PfS'] = 
i6a^d{]i.a — k^), where we have adopted the same primed notation as above. We note that 
Pa' is not necessarily 6(j)a', as it would be for a linear theory. 

4.1 Heisenberg picture and Ehrenfest's theorem 

At this point, working in the Heisenberg picture of quantum mechanics aids our discussion. 
In the Heisenberg picture, the operators are time dependent, while the states are independent 
of time. Operators follow evolution equations of the form 

- = (4.1) 

where A is some operator which has no explicit time dependence, and H is the Hamiltonian 
operator. As is well known, the explicit form of this equation for a given A is identical to the 
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evolution equation for the classical variable A, but written in terms of operators. This is true 
only up to operator ordering, and a specific ordering for the Hamiltonian must be chosen. 
Here we assume that the Hamiltonian is fully symmetric. This is known as Weyl ordering. 
When we turn to correlation functions, we will see that if we take these to be Weyl ordered 
as well, then we only ever have to deal with real valued quantities. This is a choice which can 
be made for convenience. There is no loss of useful information if this choice is made since, 
ultimately, we evaluate the statistics on super-horizon regimes where they become classical 
and the ordering is irrelevant. 
In this picture 

dx ' 1 

= Ua'pixpi + --Ua'/j'y - , (4.2) 

where Xa are the operators associated with (5i;^(q,+i)/2 s-^d 5</>q/2 for odd and even a respec- 
tively, and Ua'i3'-y' is symmetric in the last two indices. Now we wish to form evolution 
equations for the quantum correlation functions of the form 

^a'P' = '^{Xa'Xi3') , (4.3) 

where for convenience we have used the same notation as for the classical object, but now it 
is understood that Sq,/^/ is not in general symmetric in its two indices. We also have 

aa'l3'Y = {Xa'Xi^'Xy) , (4.4) 

where we have used the same notation as for the classical three-point function but now there 
is no longer symmetry in the ordering of indices. The evolution of any expectation value 
follows from Ehrenfest's theorem. In the Heisenberg picture, this amounts to Eq. (2.7), but 
with A taken to be a quantum operator. This leads to 

which is the quantum equivalent of Eq. (2.7). 

At this point similar reasoning to that in the previous sections can be applied. Iden- 
tifying (A) with Tia'i3' ill Eq. (4.3) and then with Oa'fS'y in Eq. (4.4), considering Eq. (4.5) 
and (4.2) and using Wick's theorem, we arrive at 

dX/Q,/^/ 

'/?' + Ui3' 

7' ^a'^y' ) 



dt 
dt 



= Ua'X'CtX'f^'y + Uis'yaa'X'j' + Uj'x'aa'f^'X' 

+ Ua'\'fi'^X'l3''^fi'Y + U|3'X'^l''^a'X'^^i'Y + UYX'fi''^a'X''^l3'fi' ■ (4.6) 



In deriving this equation we have used the chain rule in this quantum setting. For example, 
[xa'XjSi , H] = [xa' , H]xpi +Xai [xj3' ,H]. As before, since the quantum two and three-point func- 
tions are still accompanied by delta functions, the Fourier space integrals can be performed 
explicitly leading to analogous equations to Eqs. (2.11) and (2.13). Indeed the equations are 
identical to those earlier equations except that the indices are ordered in the same way as 
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in the equations (4.6) above. Moreover, formal solutions are possible. In this case, since 
ordering is again important, we have 



Sa'/3' = Ta'i'T 13/ jiT.il j' 



^a'i'^l3'j'^Yk'Oiiij'k' + ^a'i'j'^/S'k'^yi'^i'k'^j'l' 
^fS'i'j'^a'k'^yi'^k'i'^j'l' + Tyi'j'Ta'k'T fS'l'^k'i'^l'j' 



(4.7) 



where the T functions are identical to the classical case and follow identical evolution equa- 
tions. Therefore, the same formal integral solutions for the T functions hold, and moreover 
we find 



ix)'^fic>a-{koc, k^, ky)Ti^i3(^k^)Tiij-j(^ky) 



to 



+ r/3/i(^(3)^/i;>5-(^|3 5 ka, ky)Tiaj^[ka,)Ti0-^(ky) 

+ ^'yfi{ky)up_i^a{ky,kix,k^)Tiai::,[ka.)Ti0a-[kp,)] . 



(4.8) 



One can verify that direct differentiation of Eq. (4.8) gives the correct evolution equation. 
4.2 Connection to the In-In formalism 

In the In-In formalism, the integral expression for a correlation function (for example a 
three-point function with an initial value of zero) is 



(Xa'X/j'Xy) = — i / dt'{ Xa'Xj^iX^i ,'H\-at{t') ) 



(4.9) 



to 



where l-Lmt is the interaction Hamiltonian. Here the interaction picture (as opposed to the 
Heisenberg picture) is employed. Nevertheless, it is straightforward to verify that Eq. (4.9) 
leads directly to Eq. (4.8) using the expressions 



X a' X 131 Xy, Hint {t'] 



+ r. 



'y' J/' X X 



Xu','Hmt{t') 



(4.10) 



and 



[Xa',7iint{t')] = iUa'i>'jl'Xc,'Xfi' , (4.11) 

and employing Wick's theorem. To arrive at Eq. (4.10), we have used the first order part of 
the expression 

(4.12) 



X(^f — ^a'i'Xi' ~t~ {xi'Xji {x^iXj'^^ 



which is the quantum version of the result that the F functions can be identified with the 
coefficients of a 8N style Taylor expansion, see Refs. [20, 21] for a full discussion of the 
classical result. The result follows from substituting Eq. (4.12) into the right and left hand 
side of Eq. (4.2), and comparing with Eqs. (3.3) and (3.4). 
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4.3 Symmetrised correlation functions 

For the equations of motion for the correlation functions themselves, a simplification occurs 
if we choose to deal with fully symmetrised (Weyl ordered) correlation functions. Labelling 

5^(a'/3') = ^a'l3' > (4-13) 

and similarly the fully symmetric ordered three-point function as a^t^if^iy, we find 



dt 



U«'A'aA'/3V + ^"'^V^A'/S'^JTV ~ 3'""'^V^A'/3'^)i'7' (4-14) 

+ («' ^ /3' ^ 7') , 

where = — iSj^/j] is the imaginary part of the two point function, whose real part is 

the Weyl ordered two point function Sy^, . A possible advantage of writing the system in 
this way is that we deal only with real numbers. It is clear that as the dynamics become 
classical and the imaginary part of the two point function tends to zero, Eqs. (4.14) tend 
to Eqs. (2.10), and hence the classical transport equations are recovered from the quantum 
ones. 

5 Single field infiation 

As a simple example of the framework we have been discussing, we now consider single field 
inflation explicitly. We follow closely the study of Seery, Malik and Lyth [30]. In that study 
the authors showed how to calculate non-Gaussianity directly from the the field equations 
using a Greens function approach. The transport approach offers an alternative approach 
which also follows from the field equations. In order to keep things simple, the non-linear part 
of the field equations are truncated at leading order in slow-roll and we employ conformal 
time, such that r — )• — oo is the infinite past and r — )• the infinite future. We choose 
the basic variables to be 5(f) and 8(j)' = d5(f)/dT. Since there are two degrees of freedom the 
subscript indices run from 1 to 2, which in an abuse of notion we label (p and p respectively 
(even though 5(j)' is not in general the canonical momentum) . In the language of the transport 
formalism we have been developing, we have 

u^p{k) = 1, 

Upp{k) = -2n, (5.1) 

and 

up,,{h,k„k,) = -aV'"(0) - ^ (^likl + kl) - ^^^^1 

upppikuk^ks) = -^ ^^^'^^'^sr'"^'^ ' ^'-'^ 

where H = (da/dr)/a. Other components of the u coefficients are zero (at least to first order 
in slow-roll). These equations constitute Eqs. (15)-(16) of [30] (see also Refs. [31] and [32]). 
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Our choice of variables is the simplest choice possible, and for example implies that Tpi = 
dTffii/dT and Tpij = dF^jj/dr, and so they both satisfy second order (in time) differential 
equations. T(j,i, for example is a solution to the linear part of the second order field equation 
for 6(j). Given that it obeys the linear field equation and has the initial conditions discussed 
in §3, namely that F^^* = and T^^^ = F^p^* = 1, we can immediately identify F^p as the 
Greens function denoted by Seery et al. [30] as Gr^- 

By solving the relevant equation of motion we find, in the de-Sitter approximation for 
which a oc (— r)~^, 

k{r' - t)cos [k{T' - t)] - (1 + fcVV)sin [fc(r^ - r)] 
r0p(/c,r,r) = ' 

where, r is the later time, and r' the earlier time (which we will soon use as an intermediate 
time in an integral solution). This expression can be shown to be in agreement with Eq. (12) 
of Ref. [30]. 

Using the formal integral solution Eq. (4.7), the three-point can be calculated. Doing 
so we find identical expressions to those of Seery et al.. The solution can broken into pieces 
for each non-linear term in the original field equation. Considering the tip,/,,/, terms one finds 



Oi4,cj,4,{ki,k2,kz) 5 / dr' \r^p{ki,T,T')up^^{ki,k2,kz)T^^^{k2,T' ,T)Y^^^{kz,T\T) 

J — oo 

+ F0p(A;2, r, T')up^^{k2, ki, k^)T.^^{ki,T, T')T,^^{k3,T' , r) 

+ T^p{k3,T,T')up^^{k3,ki,k2)T,^^{ki,T,T')E^^{k2,T,T')] , (5.4) 

where the two point function is to be approximated by the de-Sitter expression, and we have 

((50(k,r)50(k2,r')) = ^^4k,T,T')6{k + , 

^^4k,T,T') = ^e-^'^^^-^'Hl + ikT)il - ikr') . (5.5) 

For example, considering the potential term in isolation {upip^ 5 —a'^V"'{(j))), we find an 
integral of the form of Eq. (18) of Ref. [30], and performing the integration^ one finds 



a^4,^{k3,ki,k2) 3 ^ 



f^\^-Qk^ + hY.k'^k, + -{- + l + ln\h7^.\}Y,kn , (5.6) 



where * indicates that background quantities are approximated by fixing them to the value 
they have when the scales of interest crossed the horizon, i G {1, 2, 3}, kt = Y2i ki, and 7 is 
Euler's constant. 

Using Eqs. (5.2), (5.3) and (5.5), together with 

^<Ppik,r,T') = ^e-^'=(---')(l + ifcr)r', (5.7) 

the relevant integrals for the remaining parts of Up^,/, and the Uppp terms can be shown to be 
of the form of Eqs. (21) and (25) of Ref. [30]. Combining these terms one finds 

a,,,{k3, k„ k2) 5 1 ^ E ^ E - \ E • (5-8) 

^Note the the integral must be regulated in the usual way to select the correct vacuum [1]. 
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This expression is in agreement with Eq. (29) of Ref. [30] and results found using the In- 
In formahsm [1, 4]. In this section we have used the analytic solutions which follow from 
the transport equations to find the three-point function of field fluctuations in single field 
inflation. This was to show directly in this simple example that the transport framework 
matches others in the literature. We could of course also have solved directly the coupled 
ODEs for the two- and three-point functions, Eqs. (4.6), using a numerical implementation 
of the equations. This would enable us to drop the slow-roll and de-Sitter approximations, 
and is an issue we will return to in future work. 



6 Prom field perturbations to ( 

It is often the case that the variables we choose to calculate with are not those which we 
finally want the correlations of. This is the case above where we calculate the three point 
function of the field fluctuation, but ultimately require the correlation functions of C. We 
now briefly describe how this is accommodated within our framework. 

Classically, C can be written in terms of any complete set of perturbations in the general 
form ^ 

C{k) = Na'Xa' + -Na'/3' {x^'Xp' - {x^'Xp')) , (6.1) 

were A''^' = {27r)^Na{k)6{'k — k^) and Na'p' = {2ir)^Nai3ik, ka, %)<J(k — k^ — k^). We stress 
these are not the 6N coefficients used in the 5N expansion, but are simply a compact way 
in which to write how at a given time is related to the fluctuations x^' at the same time. 
In the classical superhorizon regime we flnd 

(CO = Na'Np,^a,p, , (6.2) 

and 

(CCC) = NaiNpiN^iao-ipi-y' + Na'piN-^iNT-iTiai^tTipir' ■ (6-3) 

As in the previous sections, the presence of a delta function in both the deflnitions of the 
N functions and the correlations of x'^ implies that the integrations over Fourier space can 
be performed immediately leaving an expression with no convolutions present. We note that 
we are free to perform such a conversion to C, on large scales where the kernels of the 
quantities lose any Fourier dependence. 

The method of calculation is given elsewhere [21], but here we simply note that using 
separate-universe techniques, the N functions can be calculated for any cosmology in which 
the energy density is a function of the unperturbed quantities the Xa variables are pertur- 
bations of. Denoting these unperturbed quantities Xa, so that for example if xi = (5(^i, 
Xi = (f)i, one finds 

^0. = --^^^, (6.4) 



2H dXa ' 



and 



~ ~2HdXaXp ~ ~dX^a \h) dXi) + mdJC, \2H J dN dX^ Xp ' ^'^ 

where H and H are the Hubble rate and its derivative with respect to cosmic time, and are 
to be treated as functions of the underlying variables Xa- 
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6.1 Single field inflation 

Returning to the case of single field inflation under the slow-roll assumption, we have 



H = (6.6) 



and employing the slow-roll parameters 



V ' V 



Ml,^ , (6.7) 



one finds 



1 



'2t 



= (l - ^) • (6-8) 

These expression are in agreement with those in the literature and in particular with Eq. (32) 
of Ref. [30]. Using them, the bispectrum of can be calculated using Eq. (6.3), as was done 
in Ref. [30] , once again one finds agreement with the results of the In-In formalism as one 
must. 



7 The power spectrum in multi-field models 

Having laid out the sub-horizon transport formalism, and shown how it connects to both 
the In-In formalism and the field equation approach of Seery, Malik and Lyth [30], we wish 
to make one further connection to work in the literature before concluding. There is a well 
established methodology for calculating the power spectrum in multi-field models of inflation 
which includes the sub-horizon evolution. This was originally laid out in detail by Salopek, 
Bond and Bardeen in Ref. [28], and continues to be applied in numerical studies, for example 
in recent studies utilising the Pyflation numerical package [33] (see also [34] and references 
therein for related recent numerical implementations). It reflects the simple application of 
linear quantum field theory to multi-field infiation, and as such should agree with and recover 
the transport equations for Sq/j. We now show that this is indeed the case. 

The formalism of Bardeen et al. begins with the definition for the field fluctuation 
operators in terms of creation and annihilation operators, employing variables 6 (pa and 6<j)a, 
we have 

6Mt, k) = ^ac{t, fc)ae(k) + ^l^{t, k)al{-k), (7.1) 

where Uc and at are the usual creation and annihilation operators, the matrix allows for 
coupling between the fields, and in this section only, Roman numerals label the fields and 
run from 1 to n where n is the number of fields. 

Once again abusing notation by labelling the subscript indices for our u coefficients by 
(pa and Pa, one can show that the evolution equation for the ^ab matrix is given by 



^abit, k) = Up^^^{k)^cbit, k) + Up^p^{k)i>d>{t, k) , 



(7.2) 



where we sum over c from 1 to n. 

The kernel of the quantum two point function is given by 



^PaP,{t,k) = i>ac{t,k)^Ut,k). (7.3) 

Employing these relations, together with the evolution equation (7.2), to find an evolution 
equation for T,ai3{k), and recalling that Uci)^pf^{k) = 6ab and U(i,^^^{k) = we readily find that 
we arrive back at the equations of motion for Sq,^(A;) we found in §4. This can be viewed as 
a consistency check, or an alternative derivation of the transport equations for the quantum 
two-point function. 

8 Conclusions 

In this work we have shown how the transport equations for the correlation functions of 
inflationary perturbations can be extended from purely classical super-horizon equations to 
equations which are also valid on sub-horizon scales. This provides ODEs for the quantum 
statistical properties of perturbations valid on all scales. A property of these equations is 
that they do not involve a convolution over all scales, even for non-linear statistics, unlike the 
non-linear field equation for a given Fourier mode of a perturbation. In this sense evolution 
equations for statistics are simpler than those for the perturbations themselves. An alterna- 
tive system of coupled ODEs which also exhibits this feature, and allows the calculation of 
infiationary statistics, is that of the T objects, introduced in §3 

A major aim of this paper has been to show how this formulation for calculating ob- 
servables form infiation is connected to previous work. In particular we have shown how the 
quantum transport equations reduce to the classical moment transport equations for inflation 
which have been studied in detail in previous work [15, 20, 21]. Moreover, we have shown that 
integral solutions to the transport equations are equivalent to the integral equations which 
arise in the In-In formalism, and identical to the integral Greens function solutions found by 
Seery, Malik and Lyth in the inflationary context [30]. Finally, we have shown the equiv- 
alence to the method for calculating the power spectrum in multi-field inflation developed 
by Salopeck et ai. Given the previous transport studies have established the equivalence of 
the classical transport equations to the 5N formalism, this new work shows the connections 
between all methods in common use for calculating inflationary observables, both quantum 
and classical. 

The advantage of the transport route to non-Gaussianity calculations, however, is the 
presence, and straightforward origin, of the ODEs for quantum correlation functions. Ulti- 
mately we hope that these differential equations will lead to an efficient and easy to implement 
algorithm for the numerical computation of observable statistics of perturbations produced 
by infiation, accounting for the evolution from quantum vacuum to super-horizon scales. 
Numerical studies are necessary for models with non-linear terms which cannot be well ap- 
proximated within the slow- roll framework (for example [35, 36]), even by studies which go 
beyond leading order such as [37-39], as well as for general multi-field models. While some 
numerical work in the direction has already been implemented using alternative approaches 
(for example [35, 36, 40-42]), it stops short of code for the full bispectrum in multiple field 
models. This is a pressing problem given the immense number of infiationary models, and 
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the large number of parameter choices and initial conditions possible in each. Without such 
tools it is impossible to make a judgement about the viability of a particular model in the 
light of present and future data. We will therefore present a numerical implementation of 
our scheme in future work [27]. 
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